Magnetic frustration in an iron based Cairo pentagonal lattice 
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The Fe^^ lattice in the Bi2Fe409 compound is found to materiahze the first analogue of a magnetic 
pentagonal lattice. Due to its odd number of bonds per elemental brick, this lattice, subject to first 
neighbor antiferromagnetic interactions, is prone to geometric frustration. The Bi2Fe409 magnetic 
properties have been investigated by macroscopic magnetic measurements and neutron diffraction. 
The observed non-collinear magnetic arrangement is related to the one stabilized on a perfect tiling 
as obtained from a mean field analysis with direct space magnetic configurations calculations. The 
peculiarity of this structure arises from the complex connectivity of the pentagonal lattice, a novel 
feature compared to the well-known case of triangle-based lattices. 



The pentagon, a 5-edges polygon, is an old issue in 
mathematical recreation. It forms the faces of the do- 
decahedron, one of the platonic solids whose shape is re- 
produced in biological viruses and in some metallic clus- 
ters 0. The main peculiarity of this polygon is that, 
contrary to triangles, squares or hexagons, it is impossi- 
ble to tile a plane with congruent regular pentagons, the 
tilings must involve additional shapes to fill the gaps [2], 
like in the Penrose lattice (Fig. [T]). It exists however 
several possibilities of tessellation of a plane with non- 
regular pentagons, a famous one being the Cairo tessel- 
lation whose name was given because it appears in the 
streets of Cairo and in many Islamic decorations (Fig.[T]). 




FIG. 1: (color online) Pentagonal lattices, (a) Cairo pen- 
tagonal tiling, (b) Penrose pentagonal tiling, (c) pentagonal 
lattice studied in refs. (d) dodecahedron. 

Such a lattice could attract interest in the field of ge- 
ometric frustrated magnetism. Indeed, frustration usu- 
ally arises when all pairs of magnetic interactions are 
not simultaneously satisfied in a system due to the lat- 



tice topology, resulting in unusual physical properties [3[ . 
Today, most of the studied systems are mainly based on 
triangles (or tetrahedra in three dimensions), and can 
be made perfect from both structural (regular polygons) 
and magnetic (equal first neighbor magnetic interactions) 
points of view. This is not the case with pentagonal lat- 
tices that must involve only non-regular pentagons to fill 
the space. The experimental realization of such a model 
system is thus a unique opportunity that could bring 
new interesting features in this field of magnetic frus- 
tration. In this letter we show that the Fe^+ lattice in 
Bi2Fe409 materializes the first analogue of a magnetic 
pentagonal lattice, and we interpret its peculiar non- 
collinear magnetic arrangement with respect to magnetic 
frustration. 

Bi2Fe409 is a common by-product in the synthesis of 
the multiferroic compound BiFeOa, and has been claimed 
recently to display itself multiferroic properties [4]. This 
was the initial motivation of our study: this finding re- 
quires a precise determination of the symmetries in the 
magnetic ordered phase, since ferroelectricity can only 
arise in polar space groups. Such a precise determination 
has never been carried out on single crystals of Bi2Fe409, 
and only powder neutron diffraction measurements have 
been reported [5] with no definitive conclusions due to 
the complexity of this orthorhombic structure. Actu- 
ally, there are two different sites of four iron atoms in 
Bi2Fe409: Fei occupies a tetrahedral position and Fe2 
an octahedral one. In the structure, columns of edge- 
sharing Fe2 octahedra form chains along the c axis, and 
these chains are linked together by corner-sharing Fei 
tetrahedra and Bi atoms (Fig. [2]). As a result, the lattice 
formed by the different Fe^+ magnetic atoms is quite re- 
markable. From a simple examination of the structure, 
five main magnetic superexchange interactions, Ji to J5, 
can be identified (Fig. [2j). The structure can be simply 
viewed as layers perpendicular to the c axis: the mag- 
netic coupling of these layers along the c direction in- 
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FIG. 2: (color online) The Bi2Fe409 compound (space group Pbam; a = 7.965(4), b = 8.448(5), c = 6.007(3) A at 300 K). 
(left): View of the unit cell with the environment of the two Fe ions. Bi atoms are purple, Fe yellow and O red. Projection of 
the structure on the (a, c) plane (middle) and on the (a, b) plane (right). The atoms labels are in blue. Black lines materialize 
the magnetic exchange interactions. 



volves only Fe2 atoms in the octahedron columns and is 
achieved through two interactions, Ji between the atoms 
within the unit cell and J2 between the atoms of adjacent 
cells. Within a layer, each Fei interacts with its nearest 
neighbour Fei via J4 and with two nearest neighbour 
Fe2 pair (located above and below the mean plane) via 
J3 and J5. The projection of this layer along c forms 
a pentagonal lattice with three slightly different bonds 
per pentagon (Fig. [2j). This geometry is equivalent to a 
distorted Cairo lattice, the perfect one involving convex 
equilateral pentagons with equal- length sides, but with 
different associated angles (Fig. [1]). 

What should be expected from such a geometry and 
what is actually observed ? Up to now, no experimental 
studies and very few theoretical calculations are available 
in the literature . They only concern a very particu- 
lar pentagonal tiling, derived from the hexagonal lattice 
(see lattice c of Fig(T]). Calculations within the antiferro- 
magnetic Ising model on this lattice yields a disordered 
ground state with a finite entropy per spin [^Gi, ^81] , whereas 
in the Heisenberg model the classical ground state of the 
perfect pentagonal lattice has not yet been solved [7]. 
Single crystals of Bi2Fe409 were grown by the high 
temperature solution growth method using a flux of 
Bi203 [9]. The temperature and magnetic field depen- 
dences of the magnetization along the principal crystallo- 
graphic axes were measured using a SQUID magnetome- 
ter (Quantum Design) in the temperature range 2-380 K 
and in fields up to 5.5 T. The reported data, measured 
on the very same crystal used for the neutron diffraction 
study (2.5 X 2 X 1.5 mm^), were corrected for the demag- 
netizing field effect. The DC susceptibility between 300 K 
and 800 K was measured on a home made Faraday-type 
magnetic balance in a magnetic field of 0.7 T, using an as- 
sembly of several randomly oriented single crystals. The 
reported susceptibility was corrected from diamagnetism 
using the Pascal's constants Xdia = -198 10~^ emu/mol. 
Above room temperature, the magnetic susceptibility of 



Bi2Fe409 is found to obey a Curie- Weiss law with an ex- 
trapolated paramagnetic temperature Op ^ - 1670 K and 
an effective magnetic moment /ieff = 6.3(3) jiB per iron 
atom (see Fig. [3|). This later value is compatible with 
the 5.9 fiB expected for Fe^+ ions (7=5=5/2). When 
decreasing the temperature, an evident drop in the (a,b) 
plane susceptibility and a well pronounced increase of the 
c axis susceptibility are observed at Tat = 238(2) K, in- 
dicating a transition toward a magnetic long range order. 
This Neel temperature is in agreement with the results 
on single crystal reported in ref. [id but some 20 K lower 
that the one reported in refs. 0, Hi U for sintered sam- 



ples. The observed ratio between the Neel temperature 
and the paramagnetic one, Op/T^ ^ 7, suggests, as ex- 
pected from the geometry, a large degree of frustration. 

Further insights concerning the magnetic arrangement 
below Tat were obtained using neutron diffraction. The 
neutron diffraction experiments were performed on the 
two CEA-CRG diffractometers D23 and D15 at the Insti- 
tut Laue-Langevin (Grenoble, France) using wavelengths 
of respectively A = 1.280 A and 1.173A . Complete data 
collections have been made at two different temperatures 
T = 15 K and T = 300 K. The experimental data were 
corrected from absorption using the CCSL library [12|. 
The structural and magnetic arrangements were refined 
using the least-square programs MXD [13] and FULL- 
PROF [14], including extinction corrections [l5| . 

When decreasing the temperature, new reflections, 
characteristic of an antiferromagnetic order, appear 
around Tat = 244(2) K. They can be indexed with a 
propagation vector k = (1/2, 1/2, 1/2), that is a doubling 
of the magnetic cell compared to the nuclear one in the 
three directions. The temperature variation of a mag- 
netic Bragg peak is presented in Fig. [3l The magnetic 
structure was solved from the 654 magnetic Bragg peaks 
collected at T= 15 K taking into account the presence 
of two magnetic domains [9]. There are four Fei atoms 
and four Fe2 in the unit cell, whose coordinates are re- 
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FIG. 3: (color online) Linear magnetic susceptibility of 
Bi2Fe409 as a function of temperature measured in magnetic 
fields of 0.1 T and 0.7 T below and above 300 K respectively. 
Inset: Temperature variation of the (1/2, 3/2, —1/2) magnetic 
Bragg peak in zero field from neutron diffraction. 
(Bottom) Graphical representation of the magnetic structure 
refinement at T = 15 K in Bi2Fe409 : observed versus calcu- 
lated intensities. 



ported in Table HI The best refinement led to a residual 
factor Rw = 6.4%, with two nearly equally populated do- 
mains (42/58%). The magnetic moment parameters for 
the first domain are reported in Table HI The correspond- 
ing pictures of the magnetic arrangements are presented 
in Fig. |4j The moments on all the atoms are found re- 
stricted in the (a, b) plane. On each site, the four atoms 
are gathered into two pairs, antiferromagnetically cou- 
pled in the case of Fei and ferromagnetically coupled for 
Fe2. The magnetic moments of each pair are oriented 
at 90°, that is perpendicular to the moments of the 
other pair on the same site. The phase between the two 
sublattices Fei and Fe2 is 155°. The moment amplitudes 
on the two sites Fei and Fe2 are slightly different, and 
both much smaller than the expected value for a Fe^+ 
ion (5/iB)- A spin transfer process from the iron to the 
neighbouring oxygen ions [16] could explain this reduc- 
tion. 



Fei^ 0.3540 0.3361 0.5 3.52(1) 

Fei^ 0.6460 0.6639 0.5 3.52(1) 

Fei^ 0.1460 0.8361 0.5 3.52(1) 

Fei^ 0.8540 0.1639 0.5 3.52(1) 
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TABLE I: Coordinates of the Fe^+ ions in Bi2Fe409 and 
spherical components of their magnetic moments in the first 
magnetic domain. The magnetic moments are defined in 
spherical coordinates by (/x, ^,0). In cartesian coordinates 
(X,Y,Z) where X is along a, Y along b and Z along c, the mag- 
netic moment components are : Mx — fJ^ cose sine/), My — 
11 sine sincj)^ Mz = M cos(/). 



As it can be seen on Fig. IH this non-collinear mag- 
netic structure forms imbricated rectangles of spins. This 
is quite unusual in such an orthorhombic system, where 
all the directions are non-equivalent and where there is 
no fourfold axis. This peculiarity certainly has an ori- 
gin in the frustration of the magnetic interactions. One 
should also notice that for this structure the magnetic 
space group is centrosymmetric, that is non polar. This 
result is a priori not compatible with ferroelectricity, and 
this shed some doubts on the claimed ferroelectricity ob- 
served in polycrystalline Bi2Fe409 [4]. 

To check the relevance of the Cairo pentagon lattice 
model and the role of geometric frustration on the ob- 
served magnetic structure of Bi2Fe409, we have under- 
taken a two step theoretical approach based on a lo- 
calized spin model. Calculations have been made both 
for an ideal Cairo lattice (including the perfect case 
Js = J4 = J5), and in the real geometry. The angles 
between magnetic moments in a system is imposed by 
the set of magnetic interactions. In a first step, an exact 
Fourier space analysis at zero temperature is performed, 
assuming a single-k like ordering, which allows to find the 
periodicity of the magnetic structure [17]. Once the pe- 
riodicity is known, an energy minimization in real space 
gives the respective orientation of the magnetic moments 
within the magnetic unit cell. For these calculations, the 
five exchange interaction parameters Ji to J5 already de- 
scribed have been considered, using a Heisenberg Hamil- 
tonian H = j ^kSi • Sj. All the moments were 

assumed to be of equal amplitudes. 

In the Cairo pentagonal lattice, a large set of Ji values 
yield the observed (1/2, 1/2, 1/2) propagation vector. 
The doubling of the magnetic cell along c alone requires 
however J2 to be negative and Ji positive or slightly neg- 
ative. Whatever the k = (1/2,1/2,1/2) solution found. 
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FIG. 4: (color online) Magnetic structure. First and second 
magnetic domains. (Yellow): Fei sites. (Orange): Fe2 sites. 
The second domain is deduced from the first one by apply- 
ing the symmetry operation (x + 1/2,^ + 1/2,^) lost by the 
magnetic structure, taking into account the translational part 
that shifts atoms from one cell to a neighbouring one. 

the magnetic configuration is always made of antiparal- 
lel Fei pairs of moments and parallel Fe2 ones. On a 
same site, the pairs of moments are perpendicular to 
each other. This is exactly what is observed experimen- 
tally. The only parameter that could vary as a function 
of the particular values of the Ji 's is the phase angle be- 
tween the two sites Fei and Fe2. The ~ 155° observed 
value requires J3, J4, J5 to be all negative, with a ratio 
J?,/ Jb =2. 15 for I J4 1 > I J3 1 . The observed magnetic struc- 
ture of Bi2Fe409 and in particular the 90° phase angle 
between pairs of atoms at the same site, is thus recovered 
by this model with some realistic set of the five super- 
exchange interactions J^. In the Bi2Fe409 real geometry, 
the Fe^+ lattice differs from the Cairo pentagonal one by 
the presence of two stacked Fe2 atoms, instead of a sin- 
gle atom, at two vertices in each pentagon. Performing 
the calculations in the real geometry yield very similar 
results than those of the ideal case. The perpendicular 
ferromagnetic and antiferromagnetic Fei and Fe2 pairs is 
a robust recovered feature. The only difference between 
both lattices is the domain of existence of this phase that 
extends to lower IJ3I and IJ5I with respect to IJ4I in the 



real case. The Bi2Fe409 real lattice is therefore topolog- 
ically equivalent to the Cairo one after renormalization 
of the J3 and J5 interactions with respect to the J4 one. 

To summarize, we have found that the compound 
Bi2Fe409 displays a very original non-collinear magnetic 
structure, made of four Fei moments and four Fe2 ones 
forming interpenetrating patterns of four-fold spin rota- 
tions. This magnetic arrangement is very peculiar for 
a pentagonal lattice with nearest neighbour interactions 
only and is not a mere propagation of the magnetic ar- 
rangement minimizing the energy in each individual pen- 
tagon, contrary to the case of triangles based lattices. It 
actually results from two ingredients, a high geometrical 
frustration and an additional complex connectivity (two 
kinds of sites with different coordinations). This first 
materialisation of a pentagonal Cairo lattice opens new 
perspectives in the field of magnetic frustration. 
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